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Z.  IMTRODDCTZOlf 


Plates  and  shells  are  very  comnonly  used  as  structural 
members  in  construction  of  ships  sind  other  structures.  Due 
to  their  wide  usage  as  structural  members,  much  research  has 
been  performed  to  understand  their  failure  mechanisms  emd 
modes  under  different  loading  conditions.  Fracture  is  one  of 
the  most  frequent  failure  modes  in  plates  amd  shells. 

The  first  analysis  on  plates  containing  centrally 
located  through- the- thickness  cracks  and  subjected  to 
bending  loads  «ras  performed  by  Williams  [Ref.  1]  in  1951.  In 
his  analysis,  Wlllieuns  used  the  forth-order  classical  plate 
bending  theory  with  Kirchhoff  boundary  conditions.  He  showed 
that  there  are  some  dlscrepencies  between  the  stress 
distributions  for  extensional  and  bending  problems.  In 
Williaun's  bending  study,  surface  shear  stresses  on  crack 
faces  didn't  vanish.  This  physical  violation  is  attributed 
to  the  fact  that  the  Kirchoff  bending  theory  cam  not  satisfy 
the  physical  boundary  conditions  on  the  crack  surface 
conpletely. 

This  problem  was  solved  by  Knowles  amd  Wamg  [Ref .  2]  . 
They  used  Reissner/Mindlln's  sixth  order  plate  theory  which 
allows  three  boundary  conditions  at  edges.  They  solved  the 
problem  by  employing  a  singular  integral  equation.  However, 


their  solutions  were  limited  to  plates  with  very  small 
thicknesses.  Analyses  including  effects  of  the  plate 
thickness  on  the  bending  stress  distribution  were  performed 
independently  by  Hartranft  and  Sih  [Ref.  3]  and  Wang  [Ref. 
4]. 

All  the  theories  mentioned  above  assume  that  the  crack 
surface  is  stress -free  euid  straight.  In  reality,  for  a  plate 
under  a  bending  load,  there  exists  crack  closure  on  the 
compression  side  emd  crack  opening  on  the  tension  side  of 
the  crack  surface.  A  photoelastic  e^erimental  study 
performed  by  Smith  and  Smith  [Ref.  5]  showed  partial  crack 
closure.  A  significant  change  of  local  stresses  was  observed 
as  a  result  of  the  crack  closure. 

A  finite  element  analysis  was  first  undertaiken  by  Jones 
[Ref.  6]  using  Kirchhoff's  bending  theory  to  study  effects 
of  crack  closure  in  plate  bending.  He  ccmqpared  the  closure 
and  no- closure  solutions,  emd  reported  increases  in  crack 
opening  displacements  emd  Intensity  factors  due  to  existing 
crack  closure.  In  Heming's  study  [Ref.  7],  the 
Mindlin/Reissner  sixth  order  plate  theory  was  enployed, 
whj.ch  is  more  realistic  emd  allows  satisfaction  of  three 
boundary  conditions  at  the  crack  face.  These  results  were 
contradictory  to  the  results  obtained  by  Jones  [Ref.  6]  . 
Ilsming  stated  that  there  were  significemt  decreases  in  the 
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bending  stresses  and  displacements  on  the  tension  side  due 
crack  closure. 

In  those  studies,  the  line  closure  model,  which  assumed 
the  outermost  points  from  the  neutral  eucis  on  the 
conqpressive  side  of  the  crack  surface  to  be  in  contact,  was 
used.  However,  Smith  cmd  Smith's  e^eriment  [Ref.  5]  showed 
that  there  was  closure  over  cm  area  of  the  crack  surface, 
the  size  of  which  varied  depending  on  the  plate  cmd  crack 
dimensions.  A  three  dimensional  finite  element  ainalysis  %fas 
performed  by  Altirar  and  Nambisseus  [Ref.  8]  to  investigate  the 
closure  region  and  effects  on  the  stress  distribution  of 
this  closure.  A  very  crude  finite  element  mesh  was  used  in 
their  analyses. 

Finally  Kwon  [Ref.  9]  studied  the  crack  closure  effect 
for  unidirectional  con^osite  plates  using  the  line  closure 
model  and  the  Hindi in/Reissner  plate  theory.  He  found  that, 
the  crack  closure  on  the  compression  side  of  a  bent 
xinidirectional  composite  plate  increased  the  crack  opening 
displacement  and  bending  stress  on  the  tensile  side.  This 
finding  was  in  contradiction  to  that  for  an  isotropic  plate. 
Hence,  it  was  concluded  that  results  obtained  without 
considering  crack  closure  were  nonconservative  for 
unidirectional  composite  plates.  The  nonconservative 
solutions  for  unidirectional  conposite  plates  were  caused  by 
the  effect  of  treuasverse  shear. 
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In  the  present  study,  effects  of  surface  closure  on 
stresses  and  displacements  around  the  crack  tip  in  laminated 
conposite  plates  are  investigated  using  three-dimensional 
finite  element  euialyses. 
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IZ.  FINITE  ETiEMEMT  POSMULATION 


X.  OVERVIEW 

A  finite  element  program  is  developed  to  model  a  plate 
with  centrally  located  through- the- thickness  crack  and 
subjected  to  a  constcmt  bending  load.  As  shown  in  Figure  1, 
the  cartesian  coordinate  system  is  employed  in  the  model . 
Displacements  in  the  x,  y  and  z  directions  are  represented 
by  u,  V,  and  w  respectively. 


Figure  1  Geometric  representation  of  the  plate  model 
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The  plate  is  assumed  to  have  singly  supported  boundary 
conditions  along  its  edges.  Due  to  symmetries  along  the  x 
and  z  axes ,  only  a  quarter  of  the  plate  is  analyzed  to 
reduce  the  number  of  nodes  needed  in  modeling  the  plate. 
Therefore,  two  symmetric  boundary  conditions  and  two  simply 
supported  boundary  conditions  are  required  for  the  quarter 
plate.  Three-dimensional  parallel -piped  elements  are  used  in 
constructing  the  stiffness  matrix.  Bach  node  has  three 
degrees  of  freedcxn  and,  therefore,  there  are  24  dofs  per 
element  as  shown  in  Figure  2. 


Figure  2  Nodal  coordinates  of  a  typical  finite  element 
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Displacements  of  the  plate  are  functions  of  x,  y  and  z 
coordinates.  That  is, 

«=/i 

^=Mx,y,z)  (1) 

To  interpolate  the  displacements  within  each  finite  element, 
shape  functions  are  utilized.  Displacement  can  be  expressed 
as, 

*  g  8 

v=  2  N:*Ua  »=  2  N,-Va  (2) 


where  subscript  i  represents  the  node  number,  and  the 
numeric  subscripts  1,  2,  euid  3  denote  the  displacements  in 
the  X,  y,  and  z  directions,  respectively,  and  is  a 
function  of  x,  y  auid  z. 


From 

the  elasticity 

theory. 

given  below. 

du 

^  dw 

II 

^1 

„  -§SL4.§S. 

^  By  dx 

e  +  ^ 

dy 

strains  are  esqpressed  as 


(3) 
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When  written  in  the  caatrix  form. 


=  1516*24{«>24*I  (^) 

where  {e}  is  the  strain  vector  of  size  6  by  l,  {u}*  is  a 
nodal  displacement  vector  of  size  24  by  1,  and  [B]  is  a 
matrix  containing  derivatives  of  the  shape  functions. 

The  three-dimensional  generalized  Hooke's  law  is 
defined  as, 

{o}=P]{e}  (5) 

where  [D]  is  a  6  by  6  material  property  matrix,  and  {o}  and 
{e}  are  stress  and  strain  vectors.  The  stress  and  strain 

vectors  are  written  as, 

r 

{o}  =  [Ojor,  Oyy#  Xxy,  Xxzt  'Ey*] 

(6) 

{e}  =  [Exx,  8yy#  6r*,  "Yxyt  Yjcx»  Yyx]^ 

The  elastic  strain  energy  Ccui  be  written  as, 

V  V 

In  ec[uation  (7) ,  only  the  terxns  in  parenthesis  have 
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dependency  on  the  coordinate  system.  The  volume  integral  of 
these  terms  yields  the  stiffness  matrix, 

[.<]'  =/  (8) 

V 

The  equation  of  static  equilibrium  for  an  element  Cein 
be  written  as, 

=  (9) 

where  {f}^  is  the  vector  of  forces  acting  on  the  elemental 
nodes . 

B.  CQirSTITDTIVE  BOUATZON 

The  material  property  matrix  [D]  is  a  symmetric  matrix 
which  includes  elasticity  moduli  amd  Poisson's  ratios  of  the 
material.  In  the  present  work,  isotropic,  orthotropic  and 
laminated  con^osite  plates  are  examined. 

1.  ISOTROPIC  MATERIAL 

There  are  three  independent  material  constants  in 
the  [D]  matrix  for  am  isotropic  material.  Elasticity  moduli 
and  Poisson's  ratios  are  the  same  in  every  direction.  For  cm 
isotropic  material  matrix  [D]  is  given  as. 
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£•(1-0) 

£*o 

£*o 

0  0  0 

(l■^o)•(  1-2*0) 

( 1*0  )*<  1-2*0) 

(l■Hl)•(I-2•0) 

£•0 

£*(I-o) 

£*u 

HAH 

(l-*o)*{l-2*u) 

U  U  U 

£•0 

£•0 

000 

(l+o)*(  1-2*0) 

(l+o)*(l-2*o) 

(l-foHI-3*o) 

0 

0 

0 

GOO 

0 

0 

0 

0  G  0 

0 

0 

0 

0  0  G 

(10) 


2 .  ORTHOTROPIC  KRTBRIRL 

The  nxjxnber  of  Independent  material  constants 
in  the  [D]  matrix  is  six  for  an  orthoropic  material .  The 
material  properties  have  three  principal  euces.  The  material 

properties  of  the  plate  are  considered  as,  Ej/Ej^S,  20,  40, 
and  100,  G^/Ej-O.e,  Gjj/Ej-O.S,  o„«0.25,  ^md  o„»l/9.  Inverse 

of  the  [D]  matrix  is  called  the  con^liamce  matrix  amd  it  is 
given  in  [Ref.  10]  as  shown  helow. 


\/Ei  ~y>i\IE2  — Uai/^s  0  0  0 

—\)\2lE\  XIEi  ~yi32lE3  0  0  0 

— U13/J-1  ~V23lE2  \/E3  0  0  0 

0  0  0  I/G23  0  0 

0  0  0  0  I/G31  0 

0  0  0  0  0  I/G12 


where  E  denotes  inplame  elasticity  modulus,  G  denotes  shear 
modulus,  amd  numeric  sxibscripts  1,  2,  amd  3  represent 

material  principal  auces. 
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3.  LMCIKXTBD  CGMPOSITB  MXTBRIAL 

A  laminated  con^site  plate  with  [0/90],  is 
investigated.  There  are  two  [D]  matrices  in  the  model  to 
represent  the  orientations  of  0  and  90  degrees.  The  matrix 
[D]  for  the  0  degree  layer  is  the  same  as  that  for  the  90 
degree  layer  except  that  material  properties  in  the  x  emd  z 
directions  (see  Figure  1)  are  intercheuiged  for  the  90  degree 
layer.  As  a  result,  the  coaplieuice  matrix  for  the  90  degree 
layer  is  given  below, 

XIE'i  —Mx'ilEx  0  0  0 

— usj/Es  XIE^  —MxiIEx  0  0  0 

-U31/E3  -02i/£2  0  0  0 

0  0  0  I/G21  0  0  ^  ^ 

0  0  0  0  I/G13  0 

0  0  0  0  0  I/G32 

C.  STRAZN-DISPLACSNENT  SSLATZOir 

Linear  shape  functions  are  utilized  in  the  present 
finite  element  formulation.  That  is, 
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(13) 


1^4  =(f)  •©•(¥) 

•©•(!) 

Af8  =  (f)*©*(l) 

where  a,  b,  and  c  are  dimensions  of  a  finite  element  as  seen 
in  Figure  2. 

Matrix  [B] ,  which  relates  strains  to  nodal 
displacements  of  an  element  is  formed  by  taking  partial 
derivatives  of  the  shape  fxmctions.  Hence,  six  strain 
conponents  are  e3q)ressed  as  shown  below, 

8 

/=  I 
8 

/=  1 
8 

'==2  IWUfll  0^) 

/=  1 
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8  8 

e,y=  I  Z  i;lN.*Ua] 

/  =  1  1=1 

8  8 

e«=  I  |[M*i/nl+  Z 

1=1  i=l 

8  8 

Z  |[M*i/,i]+  Z  ilN>*Ua] 

/=!  i=l 

D.  STIFFSKSS  lATftZX 

The  element  stiffness  matrix  [K]  *  is  found  by  voliuoe 
integration  of  the  term  ([B]^[D][B])  as  in  the  equation  (7), 
where  x  varies  from  0  to  a,  (  x  €  [o  a] ) 

y  varies  from  0  to  b,  (  y  €  [o  b] ) 

z  varies  from  0  to  c.  (  z  €  [o  cl ) 

The  global  stiffness  matrix  is  constructed  by  adding 
element  stiffness  matrices.  In  the  present  finite  elmnent 
mesh,  the  niunber  of  nodes  in  the  x,  y,  amd  z  directions  are 
21,  9,  and  21  respectively.  It  gives  the  halfbandwidth  of 
636  and  total  degrees  of  fredom  of  11907.  A  static  pressure 
load  is  applied  at  the  nodes  ao  the  xz  plane  lying  in  the 
middle  of  the  plate. 
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111.  BOUNDARY  CONDITIONS  AND  NODSLINO  OF  CRACK  CLOSURB 


Selecting  proper  constraints  inposed  along  the  edges  of 
a  plate  is  very  important.  As  mentioned  in  [Ref.  7], 
application  of  Kirchoff's  boundary  conditions  or 
Reisner/Mindlin  boundary  conditions  results  in  contradictory 
solutions.  It  is  more  realistic  to  employ  the 
Reissner/Mindlin  theory,  which  describes  three  boundary 
conditions  along  the  edges  of  a  plate.  Therefore,  in  the 
present  models  the  Reissner/Mindlin  BCs  are  utilized. 
Because  of  the  symmetry  involved,  a  quarter  of  the  plate  is 
modeled.  The  model  has  two  symmetric  and  two  sinply 
supported  boundary  conditions.  Symmetric  BCs  are  satisfied 
by  constraining  displacements  in  the  normal  direction  to  the 
boundary  surface.  Simply  supported  boundary  conditions,  on 
the  other  hand  are  imposed  by  constraints  in  vertical 
displacements  and  tangential  displacements  at  the  boundary 
surface. 

Crack  closure  is  a  constraint  on  the  crack  face.  There 
are  three  crack  surface  models  which  are  used  in  analyses  of 
plates  with  through-the-thickness  cracks.  The  first  model 
assumes  that  the  plate  has  no  constraint  on  the  crack 
surface,  as  seen  in  Figure  3a.  This  is  called  the 
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Figure  3a  No-Closure  Model 


Figure  3c  Surface  Closure  Model 
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no- closure  model.  The  crack  surface  is  considered  as 
stress- free  €md  straigth  . 

The  assusq;)tion  of  no- closure  allows  crack  opening  to 
occur  on  the  tension  side  of  the  crack  face.  On  the  other 
side  of  the  crack  face,  deformations  of  the  same  magnitudes 
occur  in  the  opposite  direction.  Hence,  overlapping  or 
penetration  taikes  place  by  enploying  the  no- closure  model. 
However,  penetration  is  not  possible  in  the  actual  case.  In 
reality,  there  would  be  contact  over  a  region  on  the  crack 
surface.  And  this  contact  is  called  crack  closure. 

We  cam  model  the  closure  region  in  two  different  ways. 
In  the  first  model  the  outermost  nodes  on  the  compression 
side  with  respect  to  the  neutral  cucis  are  constrained  to 
avoid  penetration.  This  is  called  the  line  closure  model  and 
is  shown  in  Figure  3b.  Alternatively,  a  region  which 
contains  some  of  the  nodes  between  the  outermost  nodes  and 
neutral  aocls  on  the  compression  side  of  the  crack  face  can 
be  constrained,  as  shown  in  Figure  3c.  This  second  model  is 
called  the  surface  closure  model. 

Plate  bending  theories  permit  only  the  line  closure 
model,  because  the  theories  assume  that  plauie  sections 
remain  plain  before  and  after  bending.  The  experimental 
study  by  Smith  amd  Smith  [Ref.  5]  shows  that  crack  closure 
taJces  place  over  a  portion  of  the  crack  face.  Therefore,  the 
surface  closure  model  is  the  closest  describtion  of  the 


crack  closure.  An  iteration  process  is  necessary  to  model 
surface  closure  using  a  contact  algorithm.  Proper  number  of 
nodes  are  constrained  to  ensure  that  there  is  no  penetration 
on  the  crack  surface. 

When  using  the  finite  element  method,  it  is  possible  to 
model  the  surface  crack  closure  using  2-D  plate  bending 
elements  or  3-D  solid  elements.  If  2-D  plate  bending 
elements  are  used  in  modeling,  a  stack  of  plate  elements  are 
used  through  the  plate  thickness.  As  a  result,  the  crack 
face  can  have  a  piecewise  linear  deformation  and  the  closure 
portion  becomes  a  rectangular  shape.  With  3-D  solid 
elements,  displacements  at  the  crack  face  can  have  more 
general  deformation.  Hence,  a  general  shape  of  surface  crack 
closure  can  be  modeled.  Because  of  this  advantage,  we  use 
3-D  elements  in  the  present  study. 
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ZV.  RESULTS  AMD  DZSCUSSZOHS 


In  the  present  study,  cracked  plates  made  of  three 
kinds  of  materials  are  examined.  These  are  isotropic, 
orthotropic  and  three  layered  composite  plates  with  an 
orientation  of  [0/90/0] .  Four  different  elasticity  moduli 
ratios,  E1/S2-5,  20,  40,  and  100,  are  considered  for 

ccmqposite  aind  orthotropic  models.  Consequently,  w  amalyzed 
nine  models  with  different  material  properties.  Three  types 
of  crack  closure  modelings,  which  are  discussed  in  the 
previous  chapter,  are  ei!q>loyed  to  model  the  crack  surface. 
These  are  surface  closure,  line  closure,  and  no-closure 
models.  All  three  models  are  ccznpared  with  each  other  for 
each  cracked  plate. 

All  the  stress  euid  displacement  data  are  conqsuted  on 
the  pleuie  containing  the  crack.  A  three  dimensional 
representation  of  the  crack  surface  for  the  no- closure  model 
is  sho%m  in  Figure  4.  In  this  figure  flat  section  is  the 
xincracked  portion  of  the  crack  pleme.  Since  for  the 
no- closure  model  there  is  no  constraint  on  the  crack 
surface,  displacements  with  the  seune  magnitudes  and  in  the 
opposite  directions  are  observed  on  the  tension  emd 
compression  sides.  In  the  surface  closure  model,  due  to  the 
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igure  4  Three  Dimensional  Representation  of  the 
urface  for  No-closure  model. 
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Crack 


Figure  5  Three  Dimensional  Representation  of  the  Crack 
Surface  for  surface  closure  model. 


crack  closure,  only  displacements  in  the  crack  opening 
direction  are  allowed,  as  shotm  in  Figure  5. 

Effects  of  closure  on  the  crack  opening  displacements 
can  be  explained  by  considering  a  compressive  force  due  to 
closure  as  shown  in  Figure  6a,  that  acts  on  the  conpression 
side  of  the  crack  surface.  A  statically  equivalent  force  and 
moment,  as  in  Figure  6b,  Ceui  be  thought  acting  in  the  middle 
of  the  plane.  While  the  force  increases  the  crack  opening 
displacements  on  the  entire  surface  of  the  crack,  the  moment 
causes  decreases  on  the  tension  side.  For  smaller  elasticity 
ratios  the  effect  of  the  moment  is  greater.  Hence,  we 
observe  decreases  in  crack  opening  displacements.  For  higher 
elasticity  ratios,  especially  in  the  case  of  a  conposite 
plate,  the  effect  of  the  force  is  higher. 

When  the  elasticity  ratio  is  higher,  shear  deformation 
due  to  transverse  shear  becomes  more  inportcmt.  As  seen  in 

Figure  6c  the  rotation  caused  by  the  moment  is  decreased  due 
to  the  transverse  shear  effect.  Here  6  is  the  rotation  angle 
for  the  moment  without  considering  transverse  shear  effect, 
and  Y  is  the  rotation  angle  in  the  opposite  direction  due  to 

the  transverse  shear. 

In  the  displacement  plots,  the  horizontal  axis 
represents  the  normalized  distance  from  the  crack  tip  ouid 
the  vertical  eucis  represents  crack  opening  displacements  on 
the  tension  side  of  the  plate.  These  crack  opening 
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TENSION  SIDE 


displacements  (CODs)  are  normalized  by  the  maximum 
displacement  of  the  corresponding  no  closure  model.  Since 
there  is  contact  on  the  compression  side  for  closure  models, 
only  displacements  on  the  tension  side  are  plotted. 

It  is  observed  that  the  crack  opening  displacements  for 
line  and  surface  closure  models  are  very  close  to  each  other 
for  all  the  elasticity  ratios.  When  low  elasticity  ratios 
are  considered,  crack  closure  decreases  the  crack  opening 
displacements  on  the  tension  side  all  the  way.  This  is  due 
to  the  fact  that  the  effect  of  the  moment  is  higher  for  low 
elasticity  ratios.  Since  there  are  more  constraints  for  the 
surface  closure  model  its  displacements  are  smaller.  Figure 
7  shows  the  crack  opening  displacements  on  the  tension  side 
of  a  laminated,  cross -ply  composite  plate  with  the 
elasticity  ratio  El/E2=5. 

From  Figure  8  it  is  seen  that,  when  the  elasticity 
ratio,  E1/E2,  is  20,  crack  closure  reduces  the  displacements 
near  the  crack  tip.  On  the  other  hand,  it  increases  the 
displacements  away  from  the  crack  tip.  When  the  plots  for 
El/E2=5  and  El/E2=20  are  compared,  a  slight  increase  near 
the  crack  tip  is  also  observed. 

Further  inreases  in  the  crack  opening  displacements  are 
obtained  for  elasticity  ratio  El/E2=40,  as  seen  in  Figure  9. 
Near  the  crack  tip,  the  displacements  for  closure  cases  are 
still  lower  than  that  for  the  no  closure  case. 
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1.4 


CRACK  OPENING  OISPUCEMENTS  ON  THE  TENSION  SIDE 


Figure  7  Crack  Opening  Displacements  (CODs)  on  the  tension 
side  of  a  laminated,  cross-ply  composite  plate  with 
elasticity  ratio  E./E-=5. 
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Figure  d  Crack  Opening  Oisplacemencs  (CODs)  on  the  cension 
side  of  a  laminated,  cross-ply  coii5)osite  place  with 
elasticity  ratio  Ei/E,=20. 


crack  OPemNG  dbplacsmcnts  on  the  tension  side 


NORM.  DOT.  FROM  THE  CRACK  TIP 


Figure  9  Crack  Opening  Oisplac^ents  (CODs)  on  che  tension 
side  o£  a  laminated,  cross-ply  composite  plate  with 
elasticity  ratio  E./E,=40. 
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Figure  10  Crack  Opening  Displacements  (CODs)  on  the  tension 
side  of  a  laminated,  cross>ply  composite  plate  with 
elasticity  ratio  E,/E;=100. 


^  • 


When  E1/E2-100,  closure  increases  the  crack  opening 
displacements  all  the  way,  as  seen  in  Figure  10.  This  is 
because  the  transverse  shear  effect  is  high  enough  to 
overcome  the  effect  of  the  moment  also  near  the  crack  tip. 

Since  the  closure  cases  reduce  the  stresses  on  the 
ccmpression  side  more,  stress  distributions  on  the  tension 
side  are  different  from  those  on  the  compression  side. 
Therefore,  both  tension  and  conpression  sides  are  exaunined. 
In  the  stress  plots,  the  horizontal  aucis  represents  the 
normalized  distamce  from  the  crack  tip  upto  the  comer  of 
the  plate  along  the  uncracked  portion  of  the  crack  plame. 
The  vertical  aocis  represents  stresses  that  are  normalized  by 
the  applied  constant  pressure  . 

From  Figures  11  and  12  it  can  be  seen  that  for  B1/B2-5, 
there  exist  large  discrepencies  in  the  stress  levels  between 
closure  amd  no- closure  cases.  These  differences  occur  both 
on  the  tension  amd  the  compression  sides  of  the  plate.  On 
the  other  hamd,  the  differences  in  stresses  for  surface 
closure  amd  line  closure  cases  are  very  small  as  seen  in  the 
crack  opening  displacements.  This  is  especially  true  on  the 
tension  side. 

When  higher  elasticity  ratios  are  considered,  increases 
in  the  stresses  near  the  crack  tip  are  observed  for  closure 
cases  ccxrpared  to  no- closure  case,  as  shown  in  Figures  13  to 
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18.  When  the  E1/E2*100,  higher  stresses  are  obtained  on  the 
tension  side  for  closure  cases. 
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.  STNEiiSeS  IN  TIIK  X  OIU. 


Figure  11  Stress  distribution  on  the  tension  side  of  a 
laminated,  cross  ply  conqposite  plate  with  elasticity  ratio 
E./E.=5. 
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IN  THE  X  l)IK. 


NOKM.  STKKSSeS  IN  TIIK  X  OIR 


Figure  13  Stress  discribucion  on  the  tension  side  of  a 
laminated,  cross-ply  composite  plate  with  elasticity  ratio 
Ej/E,=20. 


n  M  'll 
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MAXIMUM  STRESSES  ON  THE  TENSION  SIDE 


NORM.  OIST.  FROM  THE  CRACX  TIP 


Figure  17  Stress  distribution  on  the  tension  side  of  a 
laminated,  cross-ply  composite  plate  with  elasticity  ratio 
E./E,=100- 


36 


Figure  18  Stress  distribution  on  the  compression  side  of  a 
laminated,  cross -ply  con^josite  plate  with  elasticity  ratio 

E1/E2-IOO. 
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V.  CONCLUSIONS  AND  RBCOHIIBNDATZONS 


When  a  laminated  composite  plate  with  a  through  the 
thickness  crack  has  a  large  ratio  between  the  longitudinal 
and  transverse  elastic  moduli,  crack  closure  increases  the 
crack  opening  displacement  and  tensile  bending  stress  around 
the  crack  tip  coitpared  to  those  obtained  without  considering 
crack  closure.  There  is  also  a  slight  difference  between 
surface  and  line  closure  solutions  since  the  surface  closure 
model  has  more  constraints.  Both  crack  closure  solutions  are 
nonconservative.  As  a  result,  crack  closure  must  be 
considered  for  most  composite  plates  "nder  consideration. 

An  experimental  study  shows  that  the  surface  closure 
model  is  the  closest  description  of  crack  closure.  Since 
similar  displacements  and  stresses  are,  however,  observed 
around  the  crack  tip,  line  closure  can  be  used  to  model 
crack  closure  for  computational  efficiency. 

In  future  studies  experiments  may  be  performed  to 
validate  the  computer  results.  Effects  of  crack  closure  on 
composite  plates  containing  delamination  cracks  need  to  be 
analyzed. 
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